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ON  TRANSONIC  FLOW  IN  HYPERBOLIC  LAVAL  NOZZLES* 
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Abstract 


Transonic  flow  in  the  throat  region  of  hyperbolic  Laval  nozzle*  is  analysed  using 
the  conform.!  curvilinear  coordinate  method.  A  solution  to  this  type  flow  is  given  for 
the  cases  of  different  radii  of  wall  curvature  at  the  throat  and  different  values  of  the 
ratio  of  the  specific  heats.  The  general  behaviour  of  the  transonic  flow  and  Mach 
number  distribution  iu  the  throat  region  are  presented.  This  method  gives  fine  accuracy. 
It  can  be  applied  in  the  designing  of  nozzles,  especially,  in  the  designing  of  aerodynamic 
laser  nozzles. 


This  paper  uses  the  conformal  curvilinear  coordinate  method 
to  analyze  the  transonic  flow  characteristics  in  the  throat  reg¬ 
ion  of  a  hyperbolic  Laval  nozzle  and  gives  a  general  solution  to 
the  transonic  flow  in  a  hyperbolic  Laval  nozzle.  This  solution 
is  applicable  to  various  radii  of  wall  curvature  at  the  throat 
and  different  ratios  of  specific  beats.  Through  calculation,  the 
relation  between  major  parameters  of  the  transonic  flow  at  the 
throat  under  various  specific  beat  ratios  and  radius  of  curvature 
as  well  as  the  typical  Mach  number  distribution  of  the  transonic 
flow.  This  method  involves  simple  calculations  and  gives  high 
accuracy.  It  can  be  used  as  a  reference  in  the  designing  of 
nozzles,  especially  aerodynamic  laser  nozzles. 
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SYMBOLS 


a  speed  of  sound.  Equation  (1) 

b  half  height  of  the  throat 
channel  of  the  nozzle. 

Equation  (1) 

f(^)  conformal  curvilinear 

coordinate  transformation 
function,  Equation  (2) 

g(C)  flow  related  function, 
Equation  (4) 

g,  throat  cross-sectional  flow 
°  related  function,  Eauation 
(9) 

H(£,n)  magnification  factor  of 
the  x  surface  relative  to 
the  c  surface.  Equations 
(2),  (7) 

M  Mach  number.  Equation  (  3) 

q  velocity.  Equation  (3) 

r  radius  of  curvature  of  throat 

wall  surface,  Equation  (7) 

u  velocity  component  in  the  C 
direction.  Equation  (3) 

v  velocity  component  in  the  n 
direction.  Equation  (3) 

x  orthogonal  linear  coordinate 
Equation  (2) 

y  orthogonal  linear  coordinate 
Equation  (2) 

z  complex  variable  of  the  orth¬ 
ogonal  linear  coordinate, 
Equation  (2) 

Y  ratio  of  specific  heats. 
Equation  (4) 

6  angle  between  the  direction  of 
the  flow  and  the  direction  £ 
Equation  (5) 


C  complex  variable  of  con¬ 
formal  curvilinear  coor¬ 
dinates,  Equation  (2) 

n  conformal  curvilinear  coor¬ 
dinate,  Equation  (2) 

nb  approach  angle  of  the  wall 

surface  of  a  hyperbolic  Laval 
nozzle.  Equation  (6) 

V  Mach  angle.  Equation  (5) 

£  conformal  curvilinear  coor¬ 
dinate,  Equation  (2) 

p  gas  density.  Equation  (3) 

i/j  flow  function.  Equation  (3) 

1 1>  max.  flux.  Equation  (10) 

s 

function.  Equation  (11) 


Subscripts : 

*  critical  parameters  when  the 
flow  velocity  equals  to  the 
local  sound  speed 

o  parameters  at  the  center 

c  parameters  at  the  points  on 

the  wall  surface  at  the  throat 
cross-section 

$;j  position  of  sound  velocity 
point 

n  partial  derivative  with 
respect  to  n 

5  partial  derivative  with 
respect  to  £ 
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In  the  design  of  nozzles,  the  dimensionless  radius  of  curva¬ 
ture  of  the  wall  surface  at  the  throat  (abbreviated  as  the  radius 
of  curvature  of  wall  surface  at  the  throat)  is  an  important  para¬ 
meter.  It  is  necessary  to  know  the  relations  between  the  maximum 
flow  of  the  nozzle,  the  shape  of  the  sonic  speed  line  and  flow 
velocity  distribution  in  the  transonic  region  with  the  radius  of 
curvature  of  the  wall  surface  at  the  throat.  But  this  problem 
has  not  been  solved.  After  selecting  the  radius  of  curvature  of 
throat  wall  surface,  the  calculation  of  the  corresponding  velocity 
distribution  in  the  transonic  region  to  be  given  as  the  initial  con¬ 
dition  for  the  calculation  of  the  supersonic  section  is  also 
another  problem  in  the  design  of  aerodynamic  laser  nozzles  which 
needs  to  be  solved  immediately. 

Presently  the  design  of  nozzles  primarily  uses  semi-empirical 
methods  which  are  reliable  to  a  certain  extent  when  the  radius  of 
curvature  of  the  throat  wall  surface  is  very  large.  But,  with  the 

reduction  of  the  radius  of  curvature  of  the  wall  surface  at  the 

throat,  the  sonic  speed  line  is  far  from  a  straight  line.  The 
assumption  of  radial  flow  pattern  in  the  throat  region  is  no  longer 
valid.  Therefore,  some  of  the  presently  available  methods,  such 
as  the  linear  sonic  speed  line  method  [1-33*  the  method  involving 
the  transformation  of  radial  flow  into  parallel  flow  [4,5] j  etc., 
are  no  longer  applicable. 

As  for  the  study  of  the  transonic  flow  in  the  throat  region 
of  the  hyperbolic  Laval  nozzle,  some  work  [6-9]  has  been  done  by 
now  which  can  be  applied  to  conditions  under  which  the  radius  of 
curvature  of  wall  surface  at  the  throat  is  larger.  [10]  used  a  time 

related  method  to  calculate  the  special  example  for  the  flow  of  a 

hyperbolic  Laval  nozzle  when  the  radius  of  curvature  of  the  wall 
surface  r  is  2.  Because  of  the  introduction  of  the  time  variable 
and  the  artificial  barrier  term,  the  relation  between  r  and  the 
transonic  flow  was  not  given. 
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Considering  that  the  transonic  flow  in  the  throat  of  the 
nozzle  is  primarily  affected  by  the  radius  of  curvature  of  the 
wall  surface  at  the  throat,  this  paper  used  the  conformal  curvi¬ 
linear  coordinate  method  to  provide  a  general  solution  to  the 
transonic  flow  in  a  hyperbolic  Laval  nozzle  in  a  unified  form. 
This  solution  is  not  only  applicable  to  various  radii  of  curva¬ 
ture  of  throat  wall  surface,  but  also  to  different  ratio  of 
specific  heats. 

I.  FLOW  FUNCTION  RELATIONS 

Let  us  consider  the  two-dimensional  steady  non-rotational 
isoentropic  flow  of  an  ideal  gas.  Let  us  use  the  half  height  of 
the  throat  b  and  the  central  parameters  a  and  p  as  the  bases  i  > 
introduce  some  dimensionless  quantities.  Based  on  the  meaning  of 
each  symbol  in  the  table  of  symbols,  a  is  added  over  the 

head  of  each  corresponding  symbol: 

x  —  x/b,  y  —  y/b,  r  —  r/b,  u  —  u/«f,  P^p/pm,  <J>  —  <P/a*P*b  (1) 


The  discussion  in  this  paper  uses  dimensionless  quantities. 
For  convenience,  the  sign  over  the  symbol  is  omitted.  Let 
x,y  be  the  orthogonal  coordinates  on  a  physical  plane  and  £,n 
be  the  curvilinear  coordinates  which  have  a  conformal  relation 
with  x,y.  Based  on  the  theory  of  the  complex  function,  there  are 
the  following  relations  between  the  coordinates  mentioned  above: 


*  —  KO,  *  —  *  +  f  ••  S  f  in,  h({,  rl)  —  |<*«/<fCI  (2) 

when  H(C,n)  is  defined  as  the  magnification  factor  of  the  z  plane 
relative  to  the  ?  plane.  Based  on  the  definition  of  the  flow 
function  <KC,n)  and  isoentropic  relation  of  the  flow,  the  velo 
city  components  u,  v  in  the  £,n  directions,  density  p  and  Mach 
number  M  have  th  following  relation  with  the  flow  function  if/  [11] 


(3) 


V  «■  -'ipf/pH,  +  K*  ' 

P<'+,)  _  X-±A  PJ  +  JL=-L  (*J  +  *»)  -  0 


Let  us  consider  that  the  streamline  of  a  compressible  flow  in 
the  nozzle  is  approximately  the  same  as  that  of  an  incompressible 
flow.  When  the  streamline  and  the  equipotential  line  of  the 
incompressible  flow  are  chosen  as  the  conformal  curvilinear  coor¬ 
dinates  £  and  n,  then  the  velocity  component  | u | > >  J v | .  From  the 
first  equation  in  (3),  we  have  I  |  >>  |  ^  |  .  Using  this  condition, 
as  an  approximation,  we  are  going  to  neglect  all  the  terms  con¬ 
taining  the  small  quantity  ip ^  in  the  flow  function  under  the  con¬ 
formal  curvilinear  coordinates  to  obtain  a  simplified  equation 
for  the  flow  function  under  the  conformal  curvilinear  coordinates. 
The  solution  to  this  simplified  flow  function  has  the  following 
integral  form  [11]: 


r  -  1  sW  I**1 
r  +  l 


djt 


where  H(£,n)  is  a  known  function  of  the  chosen  coordinates  (£,n). 
g(£)  is  a  function  to  be  determined  by  the  boundary  conditions. 
For  convenience,  g(£)  is  called  the  flow  related  function. 


In  the  curvilinear  coordinate  system  £,n  which  has  a  confor¬ 
mal  relation  with  the  orthogonal  linear  coordinates  x,y,  the  forms 
of  the  equations  of  characteristic  lines  are  invariant.  It  can 
be  written  as: 

-3^  —  «tg( 0  —  —  ft  —  *ia_‘(l/A/)  (5) 

•V 

Based  on  custom,  the  "+"  in  Equation  (.5)  represents  the  first 
family  of  characteristic  lines  and  represents  the  second 
family  of  characteristic  lines.  Thus,  from  Equation  (5)  we  can 
calculate  the  influence  line  which  is  the  second  family  character¬ 
istic  line  that  passes  through  the  wall  surface  point  at  the  cross- 
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section  of  the  throat.  The  extremum  characteristic  line  is  the 
second  family  characteristic  line  which  passes  through  the  sonic 
speed  point  on  the  center  line.  The  branch  line  is  the  first 
family  characteristic  line  which  passes  through  the  sonic  speed 
point  on  the  center  line  [11].  Notice  that  6,p  and  M  in  Equa¬ 
tion  (5)  are  all  functions  of  the  flow  function  i|>(£,n). 

II.  TRANSONIC  SOLUTION  OF  THE  NOZZLE 

The  solution  of  an  incompressible  potential  flow  in  a  hyper¬ 
bolic  nozzle  is  [12]  »  -  csc^shf ,  *  -  x  +  iy,  f  -  {  +  ir,  /  (6) 

where  the  coefficient  esc  n.  is  to  make  b 

b 

be  1.  Let  us  separate  the  real  and  the  imaginary  parts  of  Equa¬ 
tion  (6)  to  obtain  the  relations  between  the  orthogonal  linear 
coordinates  x,y  and  the  curvilinear  coordinates  Through  the 

use  of  differential  geometry  and  the  definition  of  H(£,n)  in 
Equation  (2),  we  obtain  the  expressions  for  the  radius  of  curva¬ 
ture  of  the  wall  surface  at  the  throat  and  the  magnification  fac¬ 
tor  H(C,n),  respectively,  as: 

*  —  wc  T/h  sh  5 CO  11], y  <m  esc  17*  ch  •»  ctg^ »,//({,  t\)  ■»  cic  ij^Csh1  f+cos'j;)*  (7  ) 

From  the  third  equation  in  (7),  we  can  see  that  corresponding 
to  different  nb  we  can  get  a  series  of  different  hyperbolic  nozzles 
shown  in  Figure  5.  Substituting  H(£,n)  intc.  Equation  (4),  we  get 
the  solution  of  the  flow  function  $(5>n)  of  the  hyperbolic  nozzle 
as 

«t,  i>  -  (J~L)y^T  f(g)  j;  [l  -  rteW wc.h’ew,)]*  dn  ^ 

when  £  *  0  and  n  =  Equation  (8)  gives  the  flow  ^  passing 
through  the  cross-section  of  the  throat.  At  this  time,  Is  only 
a  function  of  g(0).  In  the  supersonic  flow  in  the  nozzle,  i^b  should 
satisfy  the  condition  of  the  maximum  cross-sectional  flow  at  the 
throat,  i.e.,  (fy,/d;(0)]  —  0  [11].  Therefore,  from  Equation  (8) 

we  get 
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jt*  [1  —  ^’<0> tec'  ij]  ^1  —  y  •  -y  g’(0)»ln1)/t  Ji)  —  0 

r*  (9) 

where  g(0)  has  multiple  values.  The  minimum  real  positive  root 

is  the  flow  related  function  of  supersonic  flow  in  the  nozzle 

which  is  expressed  as  g  .  From  Equation  (8),  we  can  get  the 

s 

selected  n  by  substituting  g  into  it  and  letting  £  =  0  and  the 

s 


upper  limit  of  the  integral  n  =  ■ 

corresponding  to  the  selected  r. 


This  is  the  maximum  flow 


Based  on  the  boundary  condition  s,  when  n  =  =  'J'g-  Let 

us  bring  this  relation  into  Equation  (8)  and  obtain  the  relation 
to  determine  the  flow  related  function  g(£)  of  supersonic  flow 
in  nozzles : 

(JL~)  1  gCO  j§ *  [l  “  S  +  c°A)jr  1  -  </>,  (10  ) 


In  Equation  (10),  g(£)  has  multiple  values.  For  a  given  £, 
the  choice  of  g(£)  is  as  follows:  when  £  <  0,  g(£)  is  the  minimum 
positive  real  root  corresponding  to  the  flow  in  the  contraction 
region;  when  £  =  0,  g(£)  is  the  minimum  real  multiple  root  corres¬ 
ponding  to  the  cross-sectional  flow  at  the  throat;  when  £  >  0,  g(£) 
is  the  second  real  positive  root  corresponding  to  the  flow  in  the 
expansion  region  [11]. 

Based  on  Equations  (5),  (5),  (7 )  —  CIO )  and  by  choosing  y  -  1.^* 
this  paper  has  calculated  the  flow  characteristics  of  several  con¬ 
ditions  with  r  equal  to  10.6,  5,  2.1,  0.5  and  0.22  respectively. 

In  addition,  the  Mach  number  distributions,  influence  line,  extre¬ 
mum  characteristic  line  and  branch  line  under  the  conditions  that 
r  Is  2  and  y  is  1.667,  1.4  and  1.2  are  also  calculated.  The  results 
are  shown  in  Sections  4-9. 

Besides,  the  condition  that  y  is  2  corresponds  to  the  simulated 
gas  motion  using  the  free  surface  flow  of  water  [13].  Under  this 
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condition,  based  on  Equation  (8),  the  flow  function  iK£,n)  can  be 
integrated  into  the  following  simple  analytical  expressions: 

<Kf,  O  -  Y  vs(5)  -  y  sinVcsch  2S  ig-'Oh  §tg  ,,)g>(S)  '  ^ 

g(D  “  1 2  CSC  q*  (>/*crg  J7*)Vw(f)]  cos  y  [  a  ±  cos_,td(£)  1 
<A>)  —  2cig  ^cscli  2f  ig"'(th 

Under  such  conditions,  the  maximum  flow  ip  ,  the  cross-sec- 

s 

tional  flow  related  function  g  at  the  throat,  the  wall  surface 

s 

point  velocity  q(0,nb)  on  the  cross-section  at  the  throat  and  the 

center  point  velocity  q(0,0)  have  the  following  simple  relations 

with  r :  .  ,  ,  , 

4>,  «■  rJ(l  4* 

g,  *“  *‘(1  + 

*(0,  i/*)  —  r”«(l  +  0*(ctg_,r*)*  ^ 

*(0,  0)  —  fi(ctg~‘r 


III.  MAXIMUM  FLOW 

The  flow  discussed  in  the  section  is  the  dimensionless  flow 
expressed  by  the  critical  parameters  ?*>  »•  and  the  half  height  of 
the  throat  channel  b.  It  is  also  generally  called  the  flow  coeff¬ 
icient.  For  nozzles  with  a  given  wall  surface  pattern,  because 
of  the  limitation  of  the  cross-section  at  the  throat,  there  is  a 
maximum  flow  value  permitted  to  pass  through  the  cross-section  of 

the  throat  which  is  the  maximum  flow  i p„.  The  maximum  flow  tp  is 

s  s 

an  Important  parameter  of  nozzle  flow  and  it  is  a  necessary  condi¬ 
tion  of  supersonic  flow  in  the  nozzle.  It  is  also  the  necessary 
condition  for  the  assurance  of  the  continuous  transfer  from  sub¬ 
sonic  to  supersonic  in  the  transonic  flow  in  the  throat  region. 

In  the  study  of  the  transonic  flow  In  the  nozzle,  there  has  been 
no  good  solution  to  the  problem  of  the  maximum  flow  passing  through 
the  cross-section  at  the  throat.  For  example,  in  the  relaxation 
iteration  method,  because  the  maximum  flow  is  not  determined, 
it  was  assumed  to  be  1  [14],  This  assumption  can  only  be  applied 
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approximately  when  r  is  large.  [9]  gave  the  maximum  flow  in  the 
form  of  an  Infinite  series  which  has  poor  convergence  and  it 
showed  vibration. 

This  paper  carried  out  calculations  based  on  Equations  (8) 
and  (9).  In  Figure  1  the  curves  of  the  maximum  flow  ip  are  shown 

u 

versus  r  when  y  is  1.2,  1.4  and  2.0,  respectively.  From  the 
figure,  it  was  found  that  when  r  is  larger  and  approaching  1, 

s 

is  basically  not  affected  by  y  at  this  time.  With  decreasing  r, 
ip  becomes  smaller.  For  example,  when  y  is  1.4,  r  =  1.0,  ip  = 

S  5 

0.987;  r  =  0.50,  ip  =  0.970;  y  =  0.1,  ip  =  O.883.  For  a  given  r, 

s  s 

if  y  increases,  then  1 p  is  decreased.  The  smaller  r  is,  the 
larger  the  effect  of  y  on  ^ s . 

IV.  SONIC  SPEED  LINE 

Sonic  speed  line  is  the  boundary  line  between  the  subsonic 
region  and  supersonic  region  in  the  transonic  flow.  It  Is  also  the 
Intersection  of  the  elliptical  equation  and  the  hyperbolic  equa¬ 
tion.  Therefore,  the  determination  of  the  shape  and  position  of 
the  sonic  speed  line  has  always  been  the  major  point  in  the  study 
of  transonic  flow  in  the  nozzle. 

Based  on  Equations  (9)-(H),  Figure  2  gives  curves  showing 
the  relation  between  the  sonic  speed  point  position  y  and  r  for 
the  conditions  when  y  is  1.2,  1.4  and  2.0,  respectively,  y  de- 
creases  with  increasing  y  which  is  that  with  increasing  r  the  posi¬ 
tion  of  the  sonic  speed  point  at  the  cross-section  of  the  throat 
moves  towards  the  center  line.  When  y  increases,  then  y  also 

5 

increases  which  means  that  the  position  of  sonic  speed  point  at  the 
cross-section  of  the  throat  moves  towards  the  wall  surface. 

Figure  3  gives  the  curves  showing  the  relation  between  the 
sonic  speed  points  on  the  wall  surface  as  well  as  on  the  center 
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I 


FIGURE  1  .  Relation  between  max. 
flow  and  radius  of  curvature 


Maxrnum  flw  of  nozzle  it* 

s _ 

1.2 

1.1 

2.0 

1.0 

0.98’JO 

0.9871 

0.9811 

2.0 

0.9961 

0.9957 

0.9916 

V  r 

FIGURE  2.  Relation  between  sonic 
speed  point  with  radius  of  curva¬ 
ture  I 


Key:  1 — maximum  flow  of  the  noz2le  <J>;  2 — radius  of  curvature  of 

wall  surface  r;  3 — position  of  sonic  speed  point  at  the  cross- 
section  of  throat  y;  4 — radius  of  curvature  of  wall  surface  r 


Figure  3.  Relation  between  sonic  speed  Figure  4.  Relation  between  Mach  number 
point  and  radius  of  curvature  II  at  the  cross-section  at  the  throat  and 

radius  of  curvature 

Key:  l~position  of  sonic  speed  point  on  center  streamline  x,  ;  2— 
position  of  sonic  speed  point  on  wall  surface  x;  3 — radius  of  curva¬ 
ture  of  wall  surface  r;  *1 — Mach  no.  at  wall  surface  point  Me;  5 — 

Mach  no.  at  center  point  Mo;  6 — radius  of  curvature  of  wall  surface  r. 
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line  with  r  when  y  is  1.2,  1.4  and  2.0  respectively.  The  sonic 

speed  point  on  the  wall  surface  is  located  upstream  to  the  cross- 

section  of  the  throat,  i.e.,  x  <  0.  When  y  is  1.4,  |x  |  has  a  max- 

imum  value  of  0.267  at  r  s  1.1.  With  increasing  y,  the  sonic 

speed  point  on  the  wall  surface  moves  upstream( that  is  the  value 

of  jx  1  increases  )  .  The  sonic  speed  point  on  the  center  line  is 

located  downstream  to  the  cross-section  of  the  throat,  i.e., 

x  >  0.  When  y  is  1.4,  x  has  a  maximum  of  0.295  at  y  =  0.18. 
s  s 

With  increasing  y,  the  position  of  the  sonic  speed  point  on  the 
center  line  moves  downstream. 


Figure  5.  The  shapes 
of  the  sonic  speed  lines 


Figure  5  shows  the  shapes  and  positions  of  the  sonic  speed 

lines  under  the  conditions  that  r  is  1.4  and  y  is  10.6,  5-0,  2.0, 

1.0,  0.5  and  0.22  respectively.  From  the  figure,  it  shows  that 

when  r  >  1.1,  there  is  only  a  turning  point  on  the  center  line. 

But  when  r  <  1.1,  there  is  another  turning  point  appearing  near 

the  wall  surface.  With  decreasing  r,  the  curvature  of  the  sonic 

speed  line  becomes  larger.  Referring  to  Figure  3>  it  can  be  found 

that  the  variation  of  the  absolute  value  of  the  position  of  the 

sonic  speed  x_  on  the  wall  surface.  Starting  from  r  *  10.6,  |x  | 

gradually  becomes  larger.  When  r  *  1.1,  |x  |  has  the  maximum 

s 

value.  When  r  continues  to  decrease,  |x  (  gradually  becomes  smaller. 

s 

Another  turning  point  appears  near  the  wall  surface.  In  the  range 
where  the  sonic  speed  lines  were  plotted  in  Figure  5»  the  position 
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of  the  sonic  speed  point  on  the  center  line  x  monoclinically 

s 

decreases  with  increasing;  r. 

V.  FLOW  VELOCITY  AT  THE  CROSS-SECTION  OF  THE  THROAT 

In  the  flow  in  a  nozzle,  the  flow  velocity  distribution  at 
the  cross-section  of  t;  ■■  throat  is  mainly  affected  by  r  and  y. 

When  r  is  very  large,  the  gas  flow  velocity  can  approach  the  sonic 
speed  at  the  cross-section  of  the  throat.  But  with  decreasing  r, 
the  inhomogeneity  of  the  flow  at  the  throat  rapidly  increases. 

The  linear  sonic  speed  line  assumption  [2]  is  no  longer  applicable. 

Figure  4  shows  the  curves  showing  relations  between  the  Mach 
numbers  at  the  center  point  on  the  cross-sectional  surface  of  the 
throat  as  well  as  at  the  wall  surface  point  with  r  when  y  is  1.2, 

1.4  and  2.0,  respectively.  From  this  figure,  it  can  be  seen  that 
Mc  <  1  and  M0  >  1.  The  smaller  r  is  the  larger  the  difference 
between  M0  and  Mc  becomes.  The  increase  in  y  causes  further  inhom¬ 
ogeneity  of  the  distribution  of  Mach  number  at  the  cross-section 
of  the  throat  and  the  difference  between  M0  and  Mc  increases. 

VI.  FLOW  VELOCITY  OF  THE  CENTER  LINE 

In  order  to  explain  how  r  can  quantitatively  influence  the 
flow  velocity  distribution  on  the  center  line  of  the  nozzle.  Figure 
6  gives  the  Mach  number  distributions  for  the  conditions  that  y  is 
1.4  and  r  is  10.6,  5.0,  2.0,  1.0,  1.5  and  0.22,  respectively.  The 
figure  shows  that  when  x  is  approximately  0.37,  for  different  r  all 
the  center  line  Mach  number  distribution  curves  intersect  approx¬ 
imately  at  one  point.  The  Mach  number  M  at  this  point  is  1,08. 

This  intersection  poirtt  is  the  turning  point  of  the  center  line 
Mach  number  distribution  curves.  Near  this  turning  point,  the 
acceleration  of  the  flow  is  the  fastest.  The  distribution  of  the 
Mach  number  near  the  turning  point  is  approximately  a  straight  line. 
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Figure  6.  Center  line  Mach  number  distribution 
Key:  1 — distribution  of  center  streamline  Mach  Number  M 

The  larger  r  is,  the  approximate  straight  line  region  is  also 
larger.  When  r  decreases,  the  slope  of  the  approximately  straight 
line  increases  which  is  that  the  acceleration  of  the  flow  in¬ 
creases  . 

VII.  TRANSONIC  FLOW 

Based  on  Equations  (8)-(10),  we  obtain  the  transonic  flow  In 
the  throat  region.  The  results  are  shown  in  Figures  7-10,  When 
r  is  larger,  the  distribution  of  the  equi-Macb  number  line  in  the 
throat  region  Is  relatively  homogeneous.  With  decreasing  r  the 
curvature  of  the  Mach  number  line  becomes  more  severe  and  the  more 
inhomogeneous  the  distribution  of  the  flow  velocity  in  the  flow 
field  becomes.  The  smaller  r  Is,  the  more  concentrated  near  the 
wall  surface  at  the  throat  the  acceleration  of  the  gas  flow  becomes 
and  the  flow  velocity  near  the  center  point  becomes  slower. 

Figure  8  shows  the  equi-Mach  line  distribution  In  the  throat 
region  of  a  hyperbolic  nozzle  where  y  is  1.^  and  r  is  2,  As  a  com¬ 
parison,  in  the  same  figure,  the  results  of  Cherry  [8]  and  Serra 
[10]  are  also  plotted.  Cherry  used  a  potential  Iteration  method  to 
obtain  the  precise  solution  of  flow  in  the  nozzle  with  an  approxi¬ 
mate  hyperbolic  wall  surface  for  y  *  1,^  and  r  =  2,  and  Serra  used 
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FIGJRE  9.  Ecrui-Mach  number  line  distribution 
at  the  throat  of  the  nozzle. 


1 — hyperbolic  wall  surface:  2 — Cherry  wall  surface 


Figure  10.  Distribution  of  Equi-Mach  number  line  at  the  throat 
of  the  nozzle 

a  time  related  method  to  obtain  the  distribution  of  equi-Mach 
number  line  in  the  throat  region  of  a  hyperbolic  nozzle  when  y 
is  1.4  and  r  is  2.  From  Figure  8,  it  was  found  that  the  results 
of  sonic  speed  line  in  this  work  coincide  with  those  obtained  by 
Cherry.  For  the  equi-Mach  lines  when  M  is  0.46  and  1.90,  our 
results  agree  with  results  of  Serra.  When  M  is  0.75,  1.32  and 
1.68,  our  results  are  situated  between  the  results  of  Cherry  and 
Serra. 

VIII.  THE  EFFECT  OF  y 

In  order  to  study  the  effect  of  y  on  the  flow  in  the  nozzle, 
we  carried  out  calculations  corresponding  to  the  situation  that  y 
is  1.2,  1.4,  1.667  and  2.0,  respectively.  Figure  1  shows  the  rela¬ 
tion  between  the  maximum  flow  with  r  under  different  y.  The 

results  indicate  that  when  y  increases,  decreases.  When  r  is 

s 

larger,  the  effect  of  y  on  i|i_  is  small.  When  r  is  larger,  it  can 

s 

be  approximately  considered  that  is  not  affected  by  y.  From 

s 

Figures  2  and  3,  it  can  be  found  that  when  y  increases,  the  sonic 
speed  point  on  the  cross-sectional  surface  at  the  throat  is  moving 
towards  the  wall  surface  and  the  positions  of  the  sonic  speed  point 
on  the  wall  surface  and  the  center  line  move  toward  upstream  and 
downstream,  respectively.  That  means  the  absolute  value  of  xg  is 
increasing. 
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— r  —  1.2  - r  —  1.1  - r-1.66? 


Figure  11.  The  effect  of  ratio 
of  specific  heats  on  the  equi- 
Mach  number  line 


Figure  12.  The  effect  of  ratio 
of  specific  heats  on  the  flow  at 
the  throat 

Key:  1 — name;  2 — line  pattern; 

3 — sonic  speed  line;  4 — extremum 
characteristic  line;  5 — influence 
line;  6 — branch  line 


From  Figure  4,  it  can  be  found  that  M  increases  but  M 
decreases  with  increasing  y.  This  means  that  the  increase  in  y 
makes  the  distribution  of  Mach  number  across  the  cross-section  at 
the  throat  inhomogeneous.  When  r  is  smaller,  the  effect  of  y  is 
larger.  When  r  is  larger,  the  effect  of  y  is  smaller.  From  Fig¬ 
ure  11,  we  can  see  that,  with  increasing  y,  the  equi-Macb  number 
lines  both  upstream  and  downstream  move  towards  the  throat  region. 

In  addition,  as  the  flow  speed  approaches  the  sonic  speeds,  the  dis¬ 
placement  gradually  decreases.  The  effect  of  y  on  the  shape  and 
position  of  the  sonic  speed  line  is  small.  Figure  12  shows  the^ 
positions  of  the  sonic  lines,  influence  lines,  extremum  character¬ 
istic  lines  and  branch  lines  when  y  is  1.2  and  1.667.  The  varia¬ 
tion  of  y  has  very  small  effect  on  the  shape  and  position  of  the 
sonic  speed  lines  and  the  extremum  characteristic  lines  but  it  has 
a  large  effect  on  the  Influence  lines  and  branch  lines.  The 
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influence  region  Increases  with  increasing  y  which  is  located 
inside  the  expansion  section.  The  increase  in  this  region  is 
not  affected  by  the  streamline  pattern  of  the  wall  surface  in 
the  expansion  section. 

IX.  CONCLUSIONS 

This  paper  used  the  conformal  curvilinear  coordinate  method 
to  analyze  the  characteristics  in  the  transonic  flow  region  of 
hyperbolic  Laval  nozzles  and  gave  a  general  solution  in  a  unified 
fashion  to  the  transonic  flow  in  hyperbolic  nozzles.  This  solu¬ 
tion  is  not  only  applicable  to  various  radii  of  curvature  of  the 
wall  surface  at  the  throat  but  also  to  different  ratios  of  specific 
heats.  Through  calculations,  we  obtained  the  typical  curves  of 
the  variations  of  the  characteristics  of  the  transonic  flow  at  the 
throat  with  the  radius  of  curvature  of  the  wall  surface  at  the 
throat  under  different  ratios  of  specific  heats.  It  included  the 
maximum  flow  in  the  nozzle,  sonic  speed  line  and  the  Mach  number 
distribution  on  the  center  line,  etc.  It  also  gave  the  transonic 
flow  at  the  throat  with  typical  radius  of  curvature  of  the  wall 
surface . 

We  have  carried  out  a  comparison  between  our  result  and  those 
obtained  by  Cherry  [8]  and  Serra  [10].  The  nozzles  they  discussed 
correspond  to  the  same  situation  as  r  is  2  and  y  is  1.4  in  this 
paper.  The  results  indicated  that  there  is  good  agreement  in 
terms  of  equi-Mach  number  lines.  This  explains  that  in  addition 
to  the  capability  of  selecting  the  radius  of  curvature  of  the  wall 
surface  at  the  throat  based  on  the  need  and  of  adjusting  the  ratio 
of  specific  heats  r,  as  well  as  ease  of  calculation,  the  conformal 
curvilinear  coordinate  method  has  fairly  high  accuracy. 
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EXPERIMENTAL  TECHNIQUE  AND  EXPERIMENTAL  METHOD 


COMPRESSION  PROCESSES  OF  FREE  PISTONS  AFTER  THE  OPENING  OF 
NOZZLE  IN  A  HYPERSONIC  GUN  TUNNEL  *  1) 


Qian  Fu-xing,  Xiao  Lin-kui,  Shi  Zben-ping 
Institute  of  Mechanics,  Academia  Sinica 


ABSTRACT 


The  report  is  concerned  with  compression  processes  of  free  pistons  after  the 
openin';  of  nozzle  in  a  hypersonic  gun  tunnel.  The  nozzle  flow  will  start  at  the  time 
when  the  primary  shock  running  ahead  of  the  piston .  reaches  the  diaphrara  at  the 
entrance  of  the  nozzle  and  the  piston  starts  to  decelerate  at  the  time  when  the  shock 
reflected  from  the  nozzle  strikes  it.  The  momentum  equation  for  the  motion  of  the 
piston  and  the  continuity  equation  of  mass  of  the  gas  are  derived  respectively,  and 
the  finite-difference  calculations  are  carried  out  by  computers. 

According  to  the  results  obtained,  piston  deceleration,  overrush,  rebound,  oscilla¬ 
tions,  eqnillibrium  running  and  peak  pressure  are  discussed  for  various  given  initial  con¬ 
ditions.  It  is  shown  that  compression  processes  of  free  pistons  after  the  opening  of 
nozzle  in  a  hypersonic  gun  tunnel  are  significantly  different  from  those  for  a  closed 
nozzle.  The  results  are  compared  with  those  obtained  with  equilibrium  technique  based 
on  closed  nozzle  and  with  some  experimental  data,  and  furnish  important  information 
for  the  operation  of  hypersonic  gun  tunnels. 
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This  paper  considers  the  compression  processes  of  free  pistons 
after  the  opening  of  the  nozzle  in  a  gun  tunnel  and  establishes 
the  equation  of  motion  of  the  piston  and  equation  of  flows  in  the 
nozzle.  Furthermore,  it  uses  the  finite  difference  method  and  an 
electronic  computer  to  carry  out  a  calculation.  Based  on  the  cal= 
culated  results,  it  discusses  the  related  problems,  such  as  the 
overrush,  rebound,  oscillation,  equilibrium  running  and  peak 
pressure,  etc.,  of  the  piston  after  the  opening  of  the  nozzle 
under  various  initial  conditions.  These  results  are  compared  with 
those  obtained  with  the  conventional  equilibrium  piston  running 
technique  based  on  closed  nozzle  and  with  some  experimental 
results  and  furnished  some  meaningful  reference  data  for  the 
operation  of  gun  tunnels. 


SYMBOLS 


a.  —  speed  of  sound  in  region  i  s  - 
i  (i=l,2,3,4,5,  represent  the  ^ 
various  regions  in  Figure  1  t 
same  below)  T.  - 


-reflected  shock  wave 


-  time 


-temperature  in  region  i  (°K) 


aF  —  speed  of  sound  in  the  rear  uf  — 
region  when  the  piston  is  * 
in  equilibrium 

A  —  cross-sectional  area  of  the 

driven  section  (cm2)  u^ — 

A*  —  cross-sectional  area  of  the  u  — 
throat  passage  of  the  nozzle  P 
(cm2 

u  — 

g  —  gravitational  acceleration  P^ 

M  —  mass  of  the  piston  (kg)  x  — 

Me —  equilibrium  piston  mass  (kg) 

M  — Mach  number  of  the  main  shock., _ 

wave  xlT 

PA,PR,P  ,P.  -represent  the  gas 

pressure  in  front  of  the  pis-  p_ 
ton,  behind  the  piston,  in  y  - 
equilibrium  and  in  ^region  i, 
respectively  (kg/cm'-) 

2 

£  —  peak  pressure  (kg/cm  ) 

S  —  main  shock  wave 


--equilibrium  piston  velocity, 
also  the  flow  velocity  in  the 
stationary  chamber  in  equili¬ 
brium 

—  flow  velocity  in  region  1 

—  piston  velocity  in  uniform 
motion 

—  velocity  at  the  point  the  pis¬ 
ton  starts  to  decelerate 

--  length  coordinate  which  orig¬ 
inates  at  the  entrance  of  the 
nozzle  (Figure  1) 

■-  x  coordinate  at  the  time  the 
piston  starts  to  decelerate 

—  gas  density 

—  experimental  gas  ratio  of 
specific  heats 


I .  INTRODUCTION 


Figure  1  is  the  trajectory  of  the  piston  motion  in  an  equi- 
cross-secticn  hypersonic  gun  tunnel  and  the  x-t  wave  series  dia¬ 
grams.  When  a  diaphram  breaks,  the  high  pressure  drives  the  gas 
to  suddenly  make  the  stationary  piston  to  accelerate  and  move 
rapidly  forward  along  the  driven  section.  After  a  period  of 
accelerated  motion,  the  piston  reaches  the  maximum  velocity  and 
enters  the  uniform  motion  state.  In  the  accelerating  process,  a 
series  of  compressed  wave  is  formed  in  front  of  the  piston.  They 
further  concentrate  and  form  the  main  shock  wave  S.  The  shock 
wave  S  moves  toward  the  inlet  of  the  nozzle  at  a  velocity  higher 
than  that  of  the  piston.  When  it  reaches  the  entrance  of  the 
nozzle,  it  on  one  hand  breaks  the  secondary  diaphram  which  seals 
the  inlet  to  make  the  experimental  gas  to  begin  flow  into  the  exper 
imental  section  through  the  nozzle.  On  the  other  hand,  it  is 
simultaneously  reflected  at  the  entrance  of  the  nozzle.  The  shock 
waves  which  follow  are  repeatedly  reflected  between  the  piston  and 
the  nozzle  entrance.  When  the  first  reflected  shock  wave  Sr,  meets 
the  piston,  it  forces  the  piston  to  decelerate.  But  because  the  pi 
ton  has  a  certain  kinetic  energy,  It  still  moves  forward.  In  the 
process  of  repeated  reflection  of  shock  waves,  under  general  condi¬ 
tions  the  piston  will  "oscillate"  near  the  end  of  the  tube  which 
is  the  oscillation  motion.  It  also  produces  a  high  pressure  peak  £ 
in  the  experimental  gas.  Only  after  damping  over  a  certain  period 
of  time,  the  oscillation  will  disappear.  Afterwards,  it  reaches  a 
stable  equilibrium  state.  The  pressure  of  the  experimental  gas  is 
Pe.  The  pressure  peak  p  and  equilibrium  pressure  Pg  of  the  exper¬ 
imental  gas  as  well  as  the  corresponding  temperature  rise  are  pro¬ 
duced  due  to  the  adiabatic  compression  of  the  piston  and  the  non* 
Isoentropic  compression  of  multiple  shock  waves.  In  the  piston 
decelerating  process,  a  series  of  compressed  waves  is  produced  in 
the  driving  gas  behind  It.  They  will  concentrate  to  form  a 
reverse  shock  wave  which  moves  towards  the  driving  section. 


Figure  1.  Trajectory  of  piston  motion  and  wave  series 

1 — primary  diaphram;  2 — piston;  3 — shock  wave;  4 — second  diaphram; 

5 — nozzle 

In  the  developmental  process  of  the  hypersonic  gun  tunnel, 
piston  motion  dynamics  has  attracted  a  lot  of  research  interests 
both  theoretically  and  experimentally.  The  work  performed  in  this 
field  has  been  described  in  detail  in  [1-53.  But  all  those  studies 
were  carried  out  under  the  condition  that  the  nozzle  is  not  open. 
Actually,  as  described  above,  when  the  main  shock  wave  reaches 
the  entrance  of  the  nozzle  it  breaks  the  secondary  diaphram  and 
starts  the  flow  in  the  nozzle.  However,  these  studies  did  not  con¬ 
sider  the  effect  of  the  flow  in  the  nozzle  or  the  piston  compress- 
* 

ion  processes.  et  aL,[6,7]>  measured  the  trajectory 

of  the  piston  motion  using  the  microwave  technique  and  obtained 
the  following  under  the  conditions  with  and  without  the  flow  in 
the  nozzle.  Under  the  same  initial  conditions,  the  flow  in  the 
nozzle  has  a  relaxing  action  on  the  overrush  of  the  piston.  Knoos 
[8]  theoretically  treated  the  isoentropic  compression  processes 
which  make  the  driving  gas  in  front  of  the  piston  reach  the  peak 
pressure  in  the  driving  section  of  the  shock  tube  in  a  free  piston 
drive.  The  analytical  results  indicate  that  the  sonic  speed  bole 
at  the  primary  diaphram  has  a  relaxation  reaction  to  the  motion  of 
the  piston  under  certain  conditions.  But  the  hypersonic  gun  tunnel 
situation  is  completely  different  because  the  gas  in  front  of  the 
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pi  ston  of  the  gun  tunnel  has  the  non-isoentropic  processes  due  to 
the  presence  of  shock  waves.  In  addition,  in  the  treatment  of 
Knoos  the  pressure  of  the  driving  gas  behind  the  piston  can  be 
neglected.  It  cannot  be  neglected  for  gun  tunnels.  It  is  espe¬ 
cially  true  when  using  the  equilibrium  piston  technique. 

This  paper  will  attempt  to  treat  this  problem  theoretically 
by  considering  the  compression  processes  of  the  free  pistons  after 
the  opening  of  the  nozzle  in  a  gun  tunnel  and  then  to  establish 
the  related  equations  and  followed  by  choosing  some  typical  con¬ 
ditions  to  carry  out  a  computational  analysis  in  order  to  explain 
the  effect  of  the  existence  of  the  flow  in  the  nozzle  on  the  com¬ 
pression  processes  of  the  piston  and  other  relented  parameters. 

In  the  meantime,  the  calculated  results  are  compared  with  some 
experimental  results  to  obtain  some  meaningful  conclusions. 

II.  COMPRESSION  PROCESSES  OP  FREE  PISTONS  AFTER  THE  OPENING  OF  THE 

PISTON  IN  A  HYPERSONIC  GUN  TUNNEL 

From  the  above  description,  it  can  be  found  that  in  a  gun 
tunnel  from  the  start  to  shut-down  the  piston  has  experienced  the 
processes  of  acceleration,  uniform  velocity,  deceleration  and 
oscillation  as  well  as  reaching  the  equilibrium  state.  From  the 
instance  that  the  primary  reflected  shock  wave  St)  meets  the  pis¬ 
ton  and  make  it  decelerate,  the  motion  of  the  piston  may  have  the 
following  two  conditions:  1)  because  the  piston  has  a  relatively 
large  kinetic  energy,  the  piston  may  overrush  to  a  distance  very 
close  to  the  entrance  of  the  nozzle  to  produce  an  extremely  high 
peak  pressure.  It  then  rebounds,  oscillates  and  finally  reaches 
the  equilibrium  state;  2)  if  the  mass  of  the  piston,  the  ratio  of 
the  initial  pressures  in  the  driving  and  driven  section  and  initial 
parameters,  such  as  A#/A  can  be  properly  selected,  the  piston  may 
not  have  overrushing  and  oscillation  and  it  may  stably  enter  the 
equilibrium  state.  In  the  equilibrium  state,  the  pressures  on 
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both  sides  of  the  piston  are  equal.  The  piston  moves  at  a  uniform 
velocity  u^,,  which  is  determined  by  the  steady  flow  in  the  nozzle 
towards  the  inlet  of  the  nozzle.  Under  this  condition,  the  peak 
pressure  is  eliminated  and  the  experimental  time  is  lengthened. 
This  is  the  equilibrium  piston  technique  in  considering  flow  in 
the  nozzle.  The  gun  tunnel  should  avoid  the  first  situation  to 
the  extent  possible  and  attempt  to  operate  under  the  second  condi¬ 
tion  . 


As  described  above,  the  experimental  gas  in  front  of  the  pis¬ 
ton  has  experienced  non-isoentropic  compression  by  multiple  shock 
waves  before  reaching  an  equilibrium  state.  Also  in  the  piston 
deceleration  process,  there  is  a  series  of  expansion  waves  orig¬ 
inating  from  in  front  of  the  piston.  These  expansion  waves  will 
surpass  the  shock  waves  reflected  from  the  piston.  Besides,  due 
to  the  deceleration  of  the  piston,  the  pressure  behind  the  piston 
is  also  varying.  All  these  make  the  analysis  become  more  com¬ 
plicated. 

The  following  is  a  consideration  of  these  factors  in  sequence. 
The  analysis  by  Evans  et  al. [9]  and  Winter  has  already  indicated 
that  after  being  compressed  by  the  third  order  shock  wave  Su  or 
the  fourth  order  shock  wave  sT)  the  entropy  of  the  experimental 
gas  can  be  neglected.  However,  in  reality  as  Stalker  [1]  has 
analyzed  that  by  neglecting  the  entropy  variation  after  Sr,  only 
very  small  error  is  created.  In  addition,  Courant  and  Friedricks 
verified  through  calculation  of  a  real  example  that  because  the 
sonic  speed  of  the  gas  behind  the  third  order  shock  wave  $»,  is 
higher,  when  the  expansion  waves  in  front  of  the  piston  surpass 
that  shock  wave,  it  is  drastically  reduced  to  near  the  Mach  wave. 
Therefore,  we  treat  the  first  and  second  order  shock  wave  com¬ 
pressions  of  the  experimental  gas  as  non-isoentropic  processes  but 
consider  the  variation  of  the  state  of  gas  after  it  (i.e.,  after 
region  (5))  approximately  as  isoentropic  processes.  Therefore, 
the  pressure  in  front  of  the  piston  can  be  written  as 


Pm  -  P,(.p/fi,y 


Cl) 


As  for  t ho  pressure  variation  behind  the  piston,  we  can  make 
this  analysis.  In  the  deceleration  process,  the  slope  of  the 
compression  wave  Cdt/Jx  —  l/(«  —  «))  produced  behind  the  piston 
becomes  stecj  or  and  finally  it  forms  an  enveloping  curve.  When 
the  deceleration  of  the  piston  is  constant,  this  enveloping  curve 
is  a  parabolic  curve  which  is  the  shock  wave  curve  formed  behind 
the  piston.  Its  starting  point  is  located  on  the  negative  char¬ 
acteristic  line  created  at  the  beginning  of  the  deceleration  of 
the  piston.  The  larger  the  deceleration,  the  closer  this  starting 
point  of  the  enveloping  curve  is  to  the  starting  point  of  deceler¬ 
ation.  In  the  piston  deceleration  process,  its  back  pressure  is 
determined  by  the  motion  of  the  piston  and  the  compression  generated 
by  it.  Therefore,  when  shock  waves  occur  in  this  series  of  waves, 
it  is  not  possible  to  find  closed  solutions  to  the  motion.  Let 
us  assume  that  at  any  instance  this  pressure  is  equal  to  the  press¬ 
ure  created  by  the  piston  moving  at  constant  velocity  (equivalent 
to  instantaneous  piston  velocity),  this  pressure  pf  can  be  obtained 


from  the  normal  shock  wave  equation: 


The  equation  of  motion  of  the  piston  is 


M  <Px 
$A  di* 


(2) 

(3) 


Let  us  introduce  these  dimensionless  parameters 


and  substitute  Equations  (1)  and  (2)  into  (3),  then  we  can  obtain 
the  equation  of  motion  of  the  piston  during  the  deceleration  process 


When  the  piston  overrushes  to  reach  the  peak  pressure  and  then 
rebounds  to  point  A  in  Figure  1  (maximum  rebounding  speeds),  the 
compression  waves  behind  the  piston  will  concentrate  into  shock 
waves.  After  passing  point  A,  the  pressure  behind  the  piston  is 
replaced  by  the  pressure  produced  by  the  unsteady  motion  of  piston 
based  on  the  equilibrium  state,  which  is 


r» 


p. 


(• 


dx_ 

it 


(5) 


Therefore,  after  passing  point  A,  the  equation  of  motion  of  the 
piston  should  be  changed  into 


The  continuity  equation  of  the  flow  passing  through  the  nozzle  is 


iQ*) 

a 


r  +  1 

F5T7 


(7) 


where 


Equation  (*0  or  Equations  (6)  and  (7)  are  the  equations  of 
motion  of  the  piston  and  the  continuity  equation  of  the  flow  in 
the  nozzle  after  the  opening  of  the  nozzle  in  a  hypersonic  gun 
tunnel. 


When  the  flow  in  the  nozzle  is  not  considered,  e  =  0.  At 
this  time  we  do  not  have  to  consider  the  continuity  equation  and 
the  p  in  Equations  (4)  and  (6)  should  be  replaced  by  1/x. 


The  above  is  derived  under  the  condition  of  equal  cross- 
section.  For  the  condition  that  the  driving  segment  and  the  driven 
segment  do  not  have  equal  cross-sections,  we  must  also  consider  the 
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effect  of' area  variation.  Th_s  affect  is  primarily  shown  in  the 
variations  of  parametrs  such  as  f, ,»/,/%,**  and  r*  .  It  should 
be  pointed  out  that  the  mass  of  the  pisi-on  has  some  effect  on 
the  acceleration  process  and  consequently  on  the  parameters  such 
xK*tK*  etc*  Work  in  this  area  has  been  done  by  several  people, 
it  is  not  going  to  be  discussed  here. 

III.  NUMERICAL  COMPUTATION 


Starting  from  the  point  K  at  which  the  piston  begins  to 
decelerate,  we  use  the  finite  difference  method  to  carry  out  a 
numerical  calculation.  For  all  the  parameters  before  point  K,  we 
use  the  numerical  values  of  those  in  the  shock  tube  related  flow 
analysis.  The  parameters  at  point  K  are  also  used  as  the  initial 
parameters  of  the  computation.  The  related  parameters  in  region  (5) 
have  already  taken  the  effect  of  the  flow  in  the  nozzle  into  con¬ 
sideration. 


When  the  initial  condition  is  F  =  0, 


C8) 


Using  the  Runge-Kutta  method,  we  change  the  second  order  differ¬ 
ential  equation  into  a  set  of  first  order  simultaneous  differential 
equations  and  use  the  integration  progressing  method  to  obtain  the 
numerical  solutions  for  various  initial  conditions. 


The  piston  begins  to  decelerate  from  point  K,  it  goes  through 
overrush  to  reach  the  peak  pressure  and  then  rebounds.  When  it 
reaches  the  maximum  rebound  velocity  point  A,  it  begins  to  decel¬ 
erate  again  and  then  starts  the  oscillation  process.  From  point  K 
to  point  A,  we  use  Equations  (M)  and  (7).  After  point  A,  we  use 
Equations  (6)  and  (7)  to  carry  out  the  calculation.  In  our  calcu¬ 
lation  we  assume  the  gas  is  an  ideal  gas  and  the  driving  and  the 
driven  gases  are  all  air  or  nitrogen.  The  ratio  of  specific  heats 
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300°K. 


is  y  =  1.4  and  the  initial  temperature  Tq  = 

IV.  RESULTS  AMD  DISCUSSIONS 

Let  us  choose  that  the  length  of  the  driven  segment  is  10  m, 
and  inner  diameter  is  10  cm.  Based  on  the  various  Initial  para¬ 
meter  combinations,  listed  in  Tables  1  and  2,  we  carried  out  the 
computat ion . 


TABLE  1.  Initial  parameters  of  open-ended  nozzle 


*  naan  m  (kg> 

0.05 

0.10 

0.20 

0.30 

0.50 

0.50 

0.70  1.0 

1.5 

2.0 

Z  A* /A 

0.05 

0.10 

0.15 

0.20 

0.25 

0.35 

0.50 

J  iStiS  -YMa  M, 

2.0 

2.5 

2.5 

2.1 

• 

1 

f  />.( kg/cm1) 

300 

500 

' 

TABLE  2.  Initial  parameters  of  closed  nozzle  (A*/A=0) 


*  M(kg) 

0.05 

0.075 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.50 

0.55 

0.50 

0-60 

0.70 

1.0 

1.5 

2.0 

*  srasE'um  Mi 

2.0 

2.4 

2.5 

2.0 

7  P.(kg/c,n*) 

300 

' 

Key:  1 — piston  mass  M  (kg);  2 — cross-section  ratio  of  throat 
channel  of  the  nozzle  A*/A;  3 — Mach  number  of  the  main  shock  wave 
Mg;  4 — driving  pressure  P^(Kg/cmz);  5 — mass  of  the  piston  M  (kg); 

6 — Mach  number  of  the  main  shock  wave  M0 ;  7 --driving  pressure  Ph 
(kg/cm2)  s  4 

The  computational  results  are  plotted  in  Figures  2-4.  Figure 
2  is  the  typical  situation  of  the  variations  of  the  dimensionless 
piston  trajectory  x  and  the  dimensionless  pressure  p  with  the 
dimensionless  time  t  after  the  piston  begins  to  decelerate  from 
point  K  for  various  ratios  of  the  nozzle  throat  channel  A*/A  and 
piston  masses  M.  For  their  Mg  and  values,  the  results  qualita¬ 
tively  have  the  same  trends.  Figure  3  shows  the  relations  between 

A 

the  ratio  of  the  peak  pressure  and  the  equilibrium  pressure  P/Pe 

and  A*/A  as  well  as  M  for  various  M  values  under  the  conditions 

s 

of  both  open  and  closed  nozzles.  Figure  4  shows  the  effect  of 
throat  channel  ratio  on  the  equilibrium  piston  mass. 


m 


Fig.  2  Cont'd.  The  typical  situation  of  the  variation  of  the 
of  the  dimensionless  piston  trajectory  x  and  the  dimensionless 
pressure  p  with  the  dimensionless  t  after  the  piston  begins  to 
decelerate  from  point  K  under  various  throat  channel  ratios  of  the 
nozzle  A*/A  and  piston  masses  M. 

Key:  1  -  piston  mass  M(kg);  2  -  ration  of  throat  channel. 
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The  calculated  results,  after  considering  the  compression 
processes  after  the  opening  of  the  nozzle  and  comparing  with  those 
obtained  using  the  conventional  analysis  of  compression  processes 
of  the  piston  on  the  basis  of  closed  nozzles,  indicate  the  follow¬ 
ing  characteristics: 

1.  From  the  variation  of  x  with  time  shown  in  Figure  2,  it 
was  found  that  the  existence  of  flow  in  the  nozzle  has  a  relaxa¬ 
tion  effect  on  the  rebound  of  the  piston.  This  coincides  with 
the  results  obtained  using  the  microwave  technique  in  [7]. 

2.  From  Figure  3,  it  was  found  that,  regardless  whether  the 
nozzle  is  open  or  closed,  the  ratio  of  the  peak  pressure  and  the 
equilibrium  pressure  P/Pe  increases  with  increasing  piston  mass  M 
for  a  given  Mach  number  M^ .  For  the  open  nozzle  condition,  the 
larger  A*/A  is  the  more  serious  this  situation  becomes.  Under 

the  condition  that  there  is  a  flow  in  the  nozzle,  due  to  the  effect 
of  gas  flow  in  the  stationary  chamber  on  the  piston  motion  on  one 
hand  and  the  main  shock  wave  at  the  entrance  of  the  nozzle  cannot 
be  totally  reflected  causing  the  weakening  of  the  reflected  shock 
wave  Sfl  on  the  other  hand,  the  deceleration  of  the  piston  at  point 
K  (i.e.,  the  loss  in  kinetic  energy)  becomes  smaller.  Thus  the 
opening  of  the  nozzle  has  a  relaxing  effect  on  the  rebound  of  the 

A 

piston  which  corresponding  less  makes  P/Pe  become  larger, 

3.  From  Figure  it  is  clear  that  if  we  take  the  piston  mass 
at  P/P  =  1.05  as  the  equilibrium  piston  mass  and  then  plot  the 
relations  between  the  equilibrium  piston  mass  S  and  the  ratio  of 
the  throat  channel  of  the  nozzle  for  various  Mach  numbers  of  the 
initial  shock  wave,  we  can  easily  see  that  the  effect  of  the  exist¬ 
ence  of  flow  in  the  nozzle  on  the  equilibrium  piston  mass  is  appar¬ 
ent.  That  is  for  a  certain  M  ,  the  larger  AVA  is  the  smaller  M 

s  c 

becomes . 


Besides  these,  the  calculated  results  also  indicate  some  other 
regularities  in  the  compression  processes  of  the  piston: 

1.  From  Figure  2,  it  was  found  that  for  given  M  and  A*/A, 

s 

the  tendency  for  the  piston  to  oscillate  increases  as  M  Increases. 
But  the  time  of  exhaustion  of  the  experimental  gas  flow  is  indepen¬ 
dent  of  M.  Of  course,  the  longer  A*/A  or  M  is,  the  shorter  the 
exhaustion  time  is.  The  calculated  results  also  indicate  that  for 
the  same  A*/A  and  M,  the  larger  M  is  the  higher  the  tendency  of 
piston  oscillation  becomes.  For  the  same  Mg,  the  higher  is, 
the  weaker  this  tendency  of  oscillating  motion  becomes  with  A*/A 
and  M  unfixed. 


-  »  Pt-StUl/uB1 

0.9  -  N 

V  - -  />,-300k,/™* 

Figure  4.  The  effect  of  *  a7_  ^<^0 

the  throat  channel  on  the  ^ 

equilibrium  mass  of  the  *>  06 '  " 
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2.  From  Figure  3>  it  was  found  that  the  piston  mass  which 

A 

satisfies  the  P/Pg  <1.05  condition  has  a  finite  region  which 
indicates  that  the  requirement  of  piston  mass  in  the  equilibrium 
piston  technique  is  not  very  stringent.  It  is  especially  so  when 
M3  is  lower. 

3.  As  described  above,  after  the  piston  rebounded  beyond 
point  A  (.Figure  1),  the  pressure  behind  the  piston  Pg  should  be 
replaced  by  Equation  (5)  instead  of  Equation  (2).  However,  under 
the  condition  that  P/Pe  <  1.10,  the  piston  almost  did  not  rebound 
from  Figures  2  and  3.  Therefore,  the  replacement  is  not  necessary. 


Under  the • condit ion  that  P/P  >1.10,  if  the  piston  mass  is  not  too 
large,  usually  equilibrium  can  be  reached  after  one  rebound.  The 
calculated  result  for  this  condition  indicates  that  the  use  of 
Equation  (5)  beyond  point  A  does  not  make  too  much  difference  in 
the  calculated  results.  If  the  oscillation  motion  of  the  piston 
is  vigorous,  then  Equation  (5)  must  replace  Equation  (2)  to  carry 
out  the  calculation. 

V.  CONCLUSIONS 

Similar  to  other  treatment  carried  out  by  other  people,  we 
established  the  equation  of  motion  of  the  piston  and  the  flow  equa 
tion  in  the  nozzle  using  the  assumption  that  the  compression  pro¬ 
cesses  beyond  region  (5)  are  isotropic  processes.  Through  calcu¬ 
lation  and  analysis,  we  discussed  the  problems  related  to  the  over 
rush,  rebound,  oscillation,  equilibrium  running,  and  peak  pressure 
of  the  piston  after  the  opening  of  the  nozzle  under  various  condi¬ 
tions  of  gun  tunnel  parameters.  We  compared  our  results  with  thos 
obtained  using  the  equilibrium  piston  motion  technique  based  on 
closed  nozzles  and  those  from  experiments.  The  results  indicate 
that  the  effect  of  gun  tunnel  nozzle  flow  on  the  compression  pro¬ 
cesses  of  the  piston  is  obvious.  These  results  have  also  received 
some  experimental  verification. 

In  the  numerical  computation  of  this  work,  we  have  received 
tremendous  assistance  from  comrade  Ho  Ni-Tei.  We  wish  to  express 
our  thanks  here. 
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PREDICTION  OF  TURBULENT  BOUNDARY  LAYER  SEPARATION* 


Wu  Guo-chuan 


(Nanjin  Aeronautical  Institute) 


In  order  to  define  the  velocity  distribution  which  is 
closely  related  to  the  pressure  gradient,  A.  Buri  introduced 
the  velocity  distribution  shape  factor  f 


r  ^  8**  dU  /  D3*»\- 

“  U  dx  V  V  /  ^  j 

where  6**  is  the  momentum  thickness  (“j,  (J  ~  fj)  ^);  U  is  the 
velocity  in  the  main  flow;  u  is  the  velocity  in  the  boundary 
layer;  v  is  the  kinetic  viscosity;  x  is  the  coordinate  along 
the  direction  of  the  flow;  m  is  a  constant.  When  r>0,  it 
corresponds  to  an  accelerating  flow.  When  r<0,  it  corresponds 
to  a  decelerating  flow. 

In  the  gas  inlet  model  experiment  participated  by  the 
author,  the  result  of  the  measured  velocity  distribution  of 
model  4  near  the  separation  point  in  the  boundary  layer  is 
shown  in  Figure  1.  Based  on  the  measured  velocity  distribution 
In  the  boundary  layer,  the  shape  factor  was  calculated  to 
be  r  =  -0.062.  Based  on  the  result  measured  by  Nikuvale, 
separation  began  when  r  =  0.06  or  T  <  -0.06.  Let  us  use  r= 
-0.06  In  equation  (1)  and  make  8**  (U8**/v)m  -  Z',  then  we 
obtain 


Finding  the  derivative  of  equation  (2)  with  respect  to 
x,  we  get 
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Figure  1.  The  velocity  distribution  near  the  separation 
point  in  the  boundary  layer  of  gas  inlet  model  number  4. 


Buri  has  assumed  that  the  shear  stress  on  the  wall  sur¬ 
face  in  the  turbulent  boundary  layer  and  the  ratio  of  the  dis¬ 
placement  thickness  to  momentum  thickness  HC—aVa**)  are 
functions  of  the  shape  factor  f ; 


-  Mr  >/(i~)”  -  /iOW)- 


(4) 


/.CO 


(5) 


where  Pis  the  density  of  the  fluid  and  R**  —  us**/v 


Considering  equations  (4)  and  (5),  we  can  obtain  the 

rpl 

following  from  the  momentum  equation  of  the  turbulent  motion 


f  ir(0”J-W) 

it 


(6) 


where  FCO  -  0  +  m)/,(0  -  (2  +  m  +  (1  +  m)/,(r)]r  •  F(r)  is  very 
close  to  be  a  straight  line  (see  Figure  2)  and  it  can  be 
expressed  as  F(0  — *  —  KT  .  Buri  used  the  measured  results 


Figure  2.  The  function  F(f)  =  0.017  -  ^.15f 


of  J.  Nikuradse  together  with  the  assumption  that  m  =  1/4  to 
calculate  the  parameters  c  and  k.  For  a  decelerating  flow, 
he  obtained  that  C  =  0.017  and  k  - 


From  Figure  2,  it  was  found  that  when  r  ■»  —  0.06  KO  —  °-26. 
Plugging  this  value  into  equatio:(jj4nd  then  using  equation  (2), 
we  can  get 


'4(^)/jx~0.32/U 


(7) 


Compare  equation  (3)  with  (7),  we  get 

“  UU"/U'‘  -  5.33  - <J  (8) 

The  above  analysis  indicates  that  when  v>5.33  ,  the 
boundary  layer  will  not  separate.  But  when  a  <  5.33  ,  the 
boundary  layer  will  separate.  *-«5.33  then  is  the  condition 
under  which  the  boundary  layer  is  near  separation.  From  this 
we  can  qualitatively  determine  the  conditions  under  which  the 
boundary  layer  will  not  separate  and  also  the  shape  of  the 
velocity  distribution  in  the  main  flow. 


From  equation  (8),  we  know  that  in  order  to  satisfy  the 
condition  that  a  >  5.33  it  must  be  that  U">  o  .  In  other 
words,  the  necessary  condition  for  no  separation  in  the 
boundary  layer  is 

U"  >  o  ( 9 ) 

Therefore,  the  following  conclusions  can  be  reached:  1. 

If  the  velocity  variation  curve  in  the  main  flow  region  after 
the  maximum  is  concave  downward  (U"<0)  (see  Figure  3),  the 
boundary  layer  separation  must  occur.  2.  If  the  velocity 
variation  curve  in  the  main  flow  region  is  concave  upward 
(U" > 0)  ,  then  the  boundary  layer  may  separate  or  it  may  not 
separate.  When  o<</<5.33  ,  it  separates  and  when  a >5.33 

it  does  not  separate.  3.  If  the  boundary  layer  does  not 
separate,  then  the  velocity  variation  curve  in  the  main  flow 
region  must  be  concave  upward,  i.e.  u" > o  • 

These  three  conclusions  are  used  to  analyse  the  experi¬ 
mental  results  obtained  from  the  experiments  of  the  eight  gas 
inlet  model  the  author  participated  in  previously.  Figures  4 
and  5  are  the  distribution  curves  of  the  theoretical  axial  M 
number  (neglecting  flow  loss)  in  the  main  flow  region  of  the 
symmetrical  flow  surface  (the  plane  slightly  below  the  hori¬ 
zontal  dividing  plane)  in  the  gas  inlet  model  experiments.  The 
ordinate  is  the  axial  M  number  and  the  abscissa  is  the  distance 
z  away  from  the  lip. 

The  axial  M  number  distribution  curve  reached  a  maximum 
very  near  the  lip  for  model  No.  1.  After  the  maximum,  the 
curve  apparently  concave  upward.  Experiment  proved  that 
there  was  no  separation  region  in  Model  1  which  agrees  with  the 
third  conclusion  discussed  above.  The  distribution  curves  of 
the  axial  M  number  in  the  main  flow  region  more  or  less  appear 


Figure  3.  Velocity  variation  in  the  main  flow  region  with  and 
without  separation  in  the  houndary  layer. 


Figure  4.  Figure  5. 

0  no.  1  x  no.  3  *■  no.  6  *  no.  2  •  no.  it  0  no.  5 

•  no.  7  *  no.  8 

to  be  concave  downward.  Based  on  the  first  conclusion,  the 
boundary  layer  separation  must  occur.  This  was  verified 
experimentally.  Comparing  the  axial  M  number  distribution 
curves  in  the  main  flow  region  for  models  2,  5,  6,  7,  and  8, 
it  was  found  that  the  curve  of  model  2  had  a  very  short  region 


which  v;nr.  concave  downward  after  the  maximum.  It  was  Immedia¬ 
tely  followed  by  an  apparent  concave  upward  section.  Its 
separation  region  is  the  shortest.  The  curve  for  model  8  has 
a  long  section  concave  downward,  therefore,  the  separation 
region  drags  for  very  long.  The  curves  for  models  5,  6,  and  7 
are  in  the  middle.  The  first  portions  of  the  three  curves  all 
are  concave  downward  to  more  or  less  the  same  extent.  As  for 
the  concave  upward  portion  which  follows,  curve  no.  5  appears 
the  earliest  and  it  is  also  the  most  apparent.  Its  separation 
region  also  ends  the  earliest  among  the  three  models.  The 
concave  upward  section  in  the  latter  portion  of  curve  no.  7 
appears  last  and  it  is  also  the  least  obvious  one.  Therefore, 
its  separation  region  is  the  longest  among  the  three  models. 

Model  6  Is  situated  between  the  two  mentioned  above.  The 
axial  M  number  distribution  curve  in  the  main  flow  region  for 
Model  4  appears  to  be  obviously  concave  upward  in  a  relatively 
long  region  and  immediately  it  appears  to  be  concave  downward. 

The  downward  region  is  smaller.  It  is  predicted  that  the 
separation  region  begins  later  and  its  range  is  smaller.  Experi¬ 
mental  results  completely  verified  the  conclusions  of  the  above 
analyses.  The  distance  between  the  starting  point  of  the  sep¬ 
aration  region  and  the  lip  for  model  4  is  larger  than  those 
for  other  models.  It  is  4  m  while  for  other  models  it  is  0.5  - 
2.5  m.  Curve  no.  3  obviously  appears  to  be  concave  downward 
in  the  beginning  section.  Based  on  the  first  conclusion,  the 
boundary  layer  must  separate.  Experimental  result  indeed 
showed  separation. 

Table  I  lists  the  theoretical  analytical  predictions  and 
experimental  results.  It  shows  more  clearly  the  consistency 
bewteen  theoretical  and  experimental  results. 

The  above  analysis  indicates  that  based  on  the  variation 
of  velocity  in  the  main  flow  region  we  can  predict  whether 


Table  I 


aniflw 

*  tt 

tt  *  1 

Sift*  ^ 

Hi*., (cm)  2 

Hi*  ..(cm)  3 

1 

0 

0 

2 

2.J.  6,  7 

1.0 

3.2 

i 

0 

<M 

3.7 

4 

mam*  9 

4.0 

1.7 

5 

S.  «*17  10 

0.3 

9.1 

t 

5,  «W7  3H2ST&T  11  | 

1.0 

9.2 

7 

1.5 

12.3 

8 

|  2.  5,  6,  7AJ8  lj  2.4 

12.3  ‘ 

Key:  1.  Experimental  results,  2.  The  distance  from  the  starting 
point  of  the  separation  region  to  the  lip  X-^cm),  3.  The 

distance  from  the  end  point  of  the  separation  region  to  the  lip 
X^Ccm),  4.  Predicted  situation  in  the  separation  region  based  on 

theoretical  analysis,  5.  Model  number  of  gas  inlet  duct,  6.  No 
separation,  7.  Shortest  separation  region  among  Model  no.  2,  5, 

6,  7,  and  8,  8.  With  separation,  9.  Separation  region  starts 
late  and  has  a  small  range,  10.  Separation  region  ends  the 
earliest  among  model  no.  5>  6,  and  7,  11.  The  middle  one  in 
models  5,  6,  and  7,  12.  Longest  separation  region  among  models 
5,  6,  and  7,  13-  Longest  separation  region  among  models  2,  5,  6, 

7 ,  and  8 . 


boundary  layer  separation  will  occur  or  not  to  a  certain  extent. 
It  provides  basis  for  the  design  of  gas  inlet  duct  with  good 
flow  characteristic  and  small  stagnation  pressure  loss.  Or, 
it  can  predict  the  quality  of  the  gas  inlet  duct  already 
designed.  It  can  be  envisioned  that  the  velocity  variation  in 
the  main  flow  region  can  first  be  defined  during  the  design  in 
order  to  select  a  variation  which  corresponds  to  a  flow  with 
good  behavior  and  no  boundary  layer  separation.  The  dimensions 
and  the  geometric  shape  of  the  cross-section  can  then  be  deter¬ 
mined  based  on  the  velocity  variation. 


The  necessary  condition  for  no  separation  in  the  turbulent 
boundary  layer  obtained  in  this  paper  is  actually  not  related 
to  the  parameters  associated  with  the  gas  inlet  duct  itself. 


Therefore,  this  conclusion  is  not  only  applicable  to  inlet 
ducts,  it  also  applies  to  the  turbulent  flow  Inside  cor.rrer.ror 
and  cascades. 
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